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Quantum graphs have recently been introduced as model systems to study the 
spectral statistics of linear wave problems with chaotic classical limits. It is proposed 
here to generalise this approach by considering arbitrary, directed graphs with uni- 
tary transfer matrices. An exponentially increasing contribution to the form factor 
is identified when performing a diagonal summation over periodic orbit degeneracy 
classes. A special class of graphs, so-called binary graphs, is studied in more detail. 
For these, the conditions for periodic orbit pairs to be correlated (including corre- 
lations due to the unitarity of the transfer matrix) can be given explicitly. Using 
combinatorial techniques it is possible to perform the summation over correlated 
periodic orbit pair contributions to the form factor for some low-dimensional cases. 
Gradual convergence towards random matrix results is observed when increasing 
the number of vertices of the binary graphs. 
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I. INTRODUCTION 

Universality in spectral statistics has be established numerically and experimentally for a wide 
range of linear wave problems ranging from quantum systems (Bohigas et al (1984)) to acoustic 
(Ellegaard et al (1996)) and microwave cavities (Alt et al (1997, 1999)) in two and three dimensions 
as well as quantum maps (Saraceno and Voros (1994)) and quantum graphs (Kottos and Smilansky 
(1997, 1998)), see also Guhr et al (1998) for a recent review. The universality classes are accurately 
described by random matrix theory (RMT) even though ensemble average is not performed when 
considering spectra of individual wave problems. This fundamental puzzle is not understood until 
today and indicates that the RMT-limit is reached under more general conditions than assumed 
by Wigner, Mehta, Dyson and others (see e.g. Mehta (1991)) in the original derivation of RMT - 
results. 

A few basic facts are well established by now: wave systems, whose spectral statistics follow the 
RMT-result for Gaussian unitary or orthogonal ensemble (GUE or GOE) have in common that 

a) time propagation (discrete or continuous) is a linear, unitary transformation; 

b) the dynamics of the underlying classical system is chaotic; this implies in particular that the 
classical Perron-Frobenius operator has an isolated largest eigenvalue equal to one, positive 
Liapunov exponent and an exponentially increasing number of periodic orbits; 

c) there are no systematic periodic orbit length degeneracies other than those enforced by the 
symmetries of the classical dynamics and the unitarity of the wave propagation. 

The last point is kept vague deliberately and refers to systems which fulfill condition (a) and 
(b) but are known to deviate from RMT due to number theoretical periodic orbit degeneracies; 
examples are the cat map (Hannay and Berry (1980), Keating (1991 a,b)) and arithmetic billiards 
of constant negative curvature (Bogomolny et al (1997)). I will come back to this point in the next 
sections. 
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A direct consequence of (b) is the so-called Hannay-Ozorio de Almeida (HOdA) sum rule (Han- 
nay and Ozorio de Almeida (1984), Berry (1985)), which enables one to derive universality of the 
spectral two point correlation function in the long range limit. Considerable progress in under- 
standing the universality of spectral statistics for individual systems beyond the HOdA-sum rule 
has been made only recently by studying quantum graphs. In a series of papers Smilansky and 
coworkers demonstrated numerically that quantum graphs indeed obey RMT-statistics (Kottos 
and Smilansky (1997, 1999)); they were also able to calculated the full form factor, i.e., the Fourier 
transform of the spectral two point correlation function, in terms of periodic orbits for a specific set 
of graphs with 2x2 unitary transfer matrices (Schanz and Smilansky (1999)) and reproduced An- 
derson localisation from periodic orbit theory in a similar model (Schanz and Smilansky (1999a)). 
Deviations from universal statistical behaviour for a special set of graphs - so-called star-graphs - 
could be explained in leading order by Kottos and Smilansky (1999), a systematic way to calculate 
higher order corrections has been developed by Berkolaiko and Keating (1999). 

The main advantage in studying quantum graphs is that one can construct a wide variety of 
systems with exact periodic orbit trace formulae (in contrast to, for example, semiclassical periodic 
orbit trace formulae, see Gutzwiller (1990)). Discrete time propagation on a graph corresponds 
to a unitary transformation in terms of a finite dimensional matrix and periodic orbit lengths 
are build up by a finite number of rationally independent length segments. The exactness of the 
trace formula circumvents problems due to, for example, semiclassical errors present in periodic 
orbit trace formulae for general quantum systems with continuous classical limit. Semiclassical 
approximations do in general not preserve unitarity of the quantum propagation which leads to 
exponentially growing error terms in the long time limit (Keating (1994), Tanner (1999)). Periodic 
orbit length correlations beyond the classical HOdA-sum rules can furthermore be studied in graphs 
in detail without referring to approximations; such correlations are predicted to exist due to the 
presence of spectral universality (Argaman et al (1993)). 

The quantisation procedure for graphs chosen by Kottos and Smilansky (1997, 1999) implies 
certain restrictions on the topological structure of the graph. Solving a one-dimensional Schrodinger 
equation on the connections (or edges) between vertices with various boundary conditions calls for 
the possibility of backscattering; the underlying graph must therefore be undirected, i.e., the 
possibility to go from vertex i to vertex j implies that the reversed direction from j to i also exists. 

In the following I will broaden the picture by considering unitary matrices in general. I will 
identify a unitary matrix as a transfer matrix (or 'wave propagator') on a directed graph with 
exact periodic orbit trace formula. The corresponding classical system is, as for quantum graphs, 
given by the dynamics on a probabilistic network. Such a construction has a priori, and again like 
for quantum graphs, no semiclassical limit in the sense that the classical dynamics does not remain 
the same when increasing the matrix dimension (or the size of the graph). This is, however, not a 
prerequisite when looking at the conditions (a) - (c); one can indeed easily construct graphs and 
corresponding unitary transfer matrices which fulfill the conditions above. The main motivation 
in generalising the concept of quantum graphs lies in the possibility to study a much wider class 
of graphs including in particular directed graphs. This freedom will be used in section III, [V to 
consider a special set of graphs, so-called binary graphs. Unlike for quantum graphs, the unitary 
transfer matrix of a directed graph can not be written as a function of a wavenumber k in general 
and does not have a quantum spectrum. Like for quantum maps, one studies instead the statistics 
of the spectrum of eigen-phases of the unitary matrix. 

I will introduce some basic notations for graphs in section || and will define edge and vertex 
staying rates as well as periodic orbit degeneracy classes. An exponentially increasing contribution 
to the form factor is identified when performing a diagonal summation over degeneracy classes. I 
will then focus on balanced, directed (binary) graphs with unitary transfer matrices. The form 
factor can here be written in terms of a periodic orbit length degeneracy function. This functions 



will be derived explicitly for binary graphs with up to 6 vertices in section III. Exponentially 
increasing contributions to the form factor are identified; these contributions alternate in sign and 
balance each other in a delicate way to lead to an expression for the form factor close to the RMT 



- result. The periodic orbit form factor for graphs with up to 32 vertices is calculated in section IV 
by counting the periodic orbit degeneracies directly. Convergence of the periodic orbit expressions 
towards the RMT - result is observed for graphs with and without time reversal symmetry; this 
gives rise to the hope that a periodic orbit theory may indeed be able to resolve universality of 



spectral statistics in the limit of large vertex - numbers. 



II. GRAPHS AND UNITARY TRANSFER MATRICES 

A. Introduction and notation 

A directed graph (digraph) G consists of set of vertices V(G) connected by a set of edges E(G). 
An edge leading from a vertex i to a vertex j, (i,j G V(G)), will be denoted (ij) and the ordering 
of the pair is important. I will mainly deal with directed graphs here and will omit the specification 
'directed' in the following. The order of the graph is given by the number of vertices N = \V(G)\, 
and M = \E(G)\ is the number of edges. A graph can be characterised by its N x N adjacency 
matrix A(G) being defined here as 

f 1 if ij€E(G) . 
%i \ otherwise 

the vertices i, j G V(G) may be labeled from to N — 1 for convenience. A real or complex N x N 
matrix T(G) will be a called a transfer matrix of G if 

tij (I • > Oij = 0. 

A real transfer matrix T C '(G) which preserves probability, i.e. 

N-l 

t% = 1 Vi G V(G), tij G H (1) 

is called a classical transfer matrix in what follows. T c ' is the analogue of the classical trans- 
fer or Frobenius-Perron operator for dynamical systems with continuous configuration space and 
describes the discrete time evolution of an N - dimensional vertex density vector p according to 

N-l 

Pj {n + 1) = Pi( n )> ne]N. 

i=0 

A matrix element tfj corresponds thus to the transition probability going from vertex i to j. The 
classical transfer matrix has a largest eigenvalue equal to one; the graph is ergodic if there exists a 
walk or path from i to j for every vertex i andj. A graph is 'chaotic' (and thus necessarily ergodic) 
if the modulus of the second largest eigenvalue is smaller than one. This means, an initial density 
vector p(0) converges exponentially fast towards an equilibrium state p which is the eigenvector 
corresponding to the largest eigenvalue of T cl . 

A periodic orbit of period n on a graph is a walk on the graph which repeats after n steps. 
Each periodic orbit can be labeled in terms of a vertex symbol code (viv^ ■ ■ ■ v n ) = v given by the 
vertices G V(G) visited along the walk with ViVi+i £ E(G),Vi = 1 and v n v\ £ E{G). We 

will denote the set off all periodic orbits of period n as VO n {G). 

I will in the following focus on unitary transfer matrices T. The 'classical' dynamics correspond- 
ing to the 'wave propagation' on the graph described by the unitary matrix T is then given by 
the classical transfer matrix T cl with t°j = |ty| 2 . The unitarity of T ensures the probability con- 
servation, Eq. (|l|), for T c ' and the equilibrium state is the uniform density vector p=(l,l,...l). 
The complex non-zero matrix elements of T may be written as tij — rije lLij and one identifies 
with the length of an edge (ij) and rh = tfj is the classical transition probability. 

The conditions (a) - (c) in section || are fulfilled if the graph is chaotic and the phases are 
not rationally related apart from conditions enforced due to unitarity. The spectrum of T and 
the periodic orbits in the graph are furthermore related by an exact trace formula; the density of 
states for the eigenphases {0i}i=i t N of T is given as 



JV oo 

d(6, N)=^2 6(9 - Bi) = 7^ + -Re ^ TrT n e" me (2) 

i— 1 n— 1 

and the traces TrT™ can be written as sum over all periodic orbits of period n in the graph, i.e. 
TrT™ = J^vevo A v e lLv . The amplitude A v is the product over the transition rates r ViVi+1 along 
the path and L v corresponds to the total length of the periodic orbit. 

The spectral measure studied in more detail in this paper is the so-called spectral form factor 
K(t, N); it is the Fourier transformed of the two point correlation function 

4tt 2 

R 2 (x,N) = —(d(6)d{0 + 2irx/N)} 

and the average (.) is taken over the 9 - interval [0, 2n]. The form factor written in terms of 
periodic orbits has the form (see e.g. Tanner (1999)) 

^(r,iV) = l(|TrT"| 2 ) Ar = l( £ A v A v ,e^- L ^) Ar (3) 

v,v'£VO n 

with t taking on the discrete values r = n/N and further averaging over small intervals At is 
performed. Most periodic orbits of the graph will be uncorrelated and the corresponding periodic 
orbit pair contributions will vanish after performing the r - average. There are, however, correla- 
tions in the periodic orbit length spectrum which lead to systematic deviations from a zero mean; 
the most obvious one is between orbits which are related by cyclic permutation of the vertex code 
v. The sum over those pairs of orbits leads to the HOdA - sum rule and describes the linearised 
behaviour of K(t) for r — ► (Berry (1985)). One can immediately identify another class of exactly 
degenerated orbits on graphs; this is the set of periodic orbits which passes through each edge the 
same number of times but not necessarily in the same order. After defining the so-called edge 
staying rates qij of v as the number of times a give orbit v visits a certain edge (ij), i.e. 

n 

9ii(v) =£^'W I a?)e£(G) veVO n , (4) 
1=1 

one can write the length L v and the amplitude A v of an orbit v on a graph as 

ijeE(G) ij€E(G) 

Periodic orbits whose symbol string gives rise to the same edge staying rate vector q = 
({lij}ij£E(G)) coincide in length L v and amplitude A v ; these orbits will be called tope/logically 
degenerate. The set of all topologically degenerated orbits will be called a degeneracy class (Berko- 
laiko and Keating (1999)). The number of orbits in a given degeneracy class represented by the M 
dimensional edge staying rate vector q (with M, the number of edges of the graph) will be denoted 
the (periodic orbit length) degeneracy function P n (q;G), i.e. 

P n (q; G) = |{v e VO n \ 9tf (v) - gy, Vij G E(G)}\ . (5) 

The orbits related by cyclic permutation of the symbol code are obviously in the same degeneracy 
class. 

The traces of T which enter the density of states (||) can thus be rewritten as 

TrT" = PM)^e iL «, (6) 

q eIK„(G) 

and K„(G) C INq 7 represents the subset of the M-dimensional integer lattice INjJ* containing all the 
possible edge staying rate vectors q which correspond to periodic orbits of period n of the graph G. 
Determining the lattice K„(G) and thus the possible degeneracy classes as well as the degeneracy 



function is the main problem when studying periodic orbit length correlations on graphs. I will 
come back to this point in the next section. 

The form factor (||) can now be written as double sum over the edge rate vectors q 

K(n,N) = ±( ]T A q A q ,P n (q)P n (q')e^^-V)) AT . (7) 

q,q'eK„(G) 

A new type of diagonal contribution emerges when considering periodic orbit pairs sharing a 
common q - vector. The total contribution of topologically degenerate periodic orbit pairs, which 
obviously includes the original diagonal contributions in the HOdA - sum rule, is 

K top (n,N) = ± K p n(<l)~ eatn (8) 

q eK„(G) 

and at > in general; (the rate at can be calculated using large deviation techniques (Dcmbo and 
Zcitouni (1993)), strict upper and lower bounds are < at < ht, and ht is the topological entropy 
for the graph). All the contributions to K top are positive which coincides with a result obtained 
by Whitney et al (1999) using diagrammatic techniques for periodic orbit formulae. The diagonal 
approximation Kt op ~ following from the HOdA - sum rule is valid only for small r = when 
cyclic permutation is the main source of degeneracies. (This is in general the case for those n 
values for which the majority of orbits visits a given edge at most once). 

Unitarity of the underlying T matrix implies the asymptotic result linv^oo K(r, N) — 1; the 
exponentially increasing topological contributions K top must therefore be counterbalanced by ad- 
ditional correlations in the periodic orbit length spectrum. We will show that these kind of correla- 
tions originate from the unitarity of the T matrix and that the cancelation mechanism is extremely 
sensitive leaving little space for approximate or asymptotic treatments. 

All what has been said so far is true for arbitrary unitary matrices, and thus especially for 
transfer matrices of quantum graphs and also for general quantum maps. In order to study the 
phenomena of periodic orbit correlations due to unitarity more closely, I will now focus on a special 
class of chaotic graphs with uniform transition probabilities for which all relevant periodic orbit 
correlations can be given explicitly. 

B. Binary graphs and periodic orbit correlations 

One of the simplest, non-trivial class of graphs are balanced, directed binary graphs B^; these 
are connected graphs with N vertices (N even) for which each vertex has exactly two incoming 
and two outgoing edges. The adjacency matrix Ajv of a binary graph can be written in the form 

a -I §2i ' j + §2i+1 ' j for < z < f () 

and the number of edges of is M — 2N. Some examples together with their adjacency matrices 
are shown in Figs. ^, |], and ||. It will sometimes be useful to switch from a vertex code to an edge 
code. A suitable choice is to assign each edge ij corresponding to a non-zero matrix element of 
the adjacency matrix (||) an edge code 

i e = 2i + j mod2, i e = 0, 1 . . . , 2N - 1 (10) 

The edge code is given for the examples Figs. ^ [|, and |j. 

Transfer matrices of binary graphs have been studied in connection with combinatorial problems 
for binary sequences (Stanley (1999)), as well as the semiclassical quantisation of the Anisotropic 
Kepler Problem using binary symbolic dynamics (Gutzwiller (1988), Tanner and Wintgen (1995)) 
and have been discussed in the context of general quantum maps (Bogomolny (1992)). Saraceno 
(1999) recently proposed a quantisation scheme for the baker map which also leads to quantum 
maps of the form @ . 

Binary graphs with adjacency matrices (|^) are connected, i.e., each vertex can be reached from 
every other vertex, here after at least [log 2 N] + 1 steps. The topological entropy h t = log 2 



independent of the order of the graph. The subset of binary graphs of order N — 2 k ,k £ IN, 
the so-called de Bruijn - graphs (Stanley (1999)), deserves special attention; the dynamics on 
these graphs can directly be related to the set of all binary sequences and there exists a one-to-one 
relation between finite binary symbol strings (a,\, a%, . . . a n ), a* £ {0, 1} of length n and the periodic 
orbits of the graph i.e. 

(ai,a 2 , ...On) <-> {vi,v 2 , ...v n ) a, £ {0, 1}, £ {0, 2 k - 1} 

with 

k 

Vi = a l+ j-i2 k ~ and a i+n = a 4 
i'=l 

for graphs of order N — 2 k . The number of orbits of period n on these graphs is exactly 2™. 

I will consider unitary transfer matrices of binary graphs next. The unitarity conditions for a 
transfer matrix T^r of a binary graph with adjacency matrix (^) can be stated simply in terms of 
the unitarity conditions for the set of N/2 different 2x2 matrices with 

Ui=f? ,W ?' 2i+1 1 i= 1,2,..., N/2 

V ,2i *i+f ,2i+l / 

Focusing on unitary binary transfer matrices with uniform local spreading, i.e., setting \tij \ = l/v2 
for all non-zero matrix elements of TV, the unitarity condition can be written as a pure phase 
correlation. One obtains the following relation between the lengths of edges 

+ L. 1+ n i2i+1 ) - + i i+ f >2i ) mod 2^ = 7T i = 0, 2, . . . , N/2 - 1 , (12) 

the corresponding local network is shown in Fig. [l]. The unitary condition ( [12] ) will be shown to be 
responsible for the periodic orbit correlations relevant to balance out the exponentially increasing 
topological contributions to the form factor K(t,N). Its simplicity makes it possible to turn the 
problem of finding periodic orbit length correlations into a combinatorial problem of finding all 
exact periodic orbit degeneracies (up to phase differences being a multiple of 7r), which can be done 
in principle. 

2i 




• i+N/2 



2i+l 

FIG. 1. Local network of correlated edge lengths, see Eq. (^); opposite edges form a pair, the two pairs 
have a combined length difference of ir. 

The dynamics described by the classical transfer matrix with transition probabilities tfj = 1/2 
for all possible transitions ij in the binary graph is maximally mixing for the geometry (^|), i.e. 
ht = K = log 2 and K is the Kolmogorov entropy for graphs (Schuster (1989)); the conditions (a) - 
(c) in section | are thus satisfied as long as there are no systematic edge length correlations present 
except from those introduced through Eq. (|l2"|). One can furthermore show that the generalised 
diagonal contribution (|^) increases exponentially with a rate at — h t — log 2 independent of the 
order of the binary graph. 



(11) 



Periodic orbit correlations beyond topological degeneracies can be expressed in terms of edge 
and vertex staying rates. The vertex staying rates <7j(v) of an orbit v are denned analogue to (Q) 
as the number of times a periodic orbit visits a vertex i, i.e. 



l=l,n 

Vertex and edge staying rates are connected by conservation laws (or shift invariance properties 
(Dembo and Zeitouni (1993))) of the form 

N 

qi,2i+qi,2i+i=q [ i ]t i + q [ i ]+ M. ! i = qi Vi = 0, . . . , — - 1 (14) 

N 

9i,2i-f +9i,2i-f +i =Q[±],i+Q[$]+%,i=4i Vi = —,..., iV- 1 



incoming edges outgoing edges 

and [.] denotes the integer part. A direct consequence of (|l^) and (R) is the following condition 
for periodic orbit correlations: 

All periodic orbits having the same vertex staying rates q = (<?o, ■ ■ ■ j Qn—i) differ in length exactly 
by a multiple of it . 

This can be shown by noting that for two orbits v, v' € VO n (B]y) with Aq = q(v) — q(v') = 0, 
one obtains 

A qi ,2i + Aq^i+i = 0, Aq it2 i + Aq i+ M. t 2i = °> A ^+f )2 i+i + A Qi+%,2i = > 
see also Fig. ||. One therefore has 

Ag ii2 i = Aq l+ K.2i+i = -Agj i2 i+i = -Ag^^ 
which together with (jlj) yields 

N/2-1 

AL = L v — L v / = 7r Agr ij2 i = 7rd ViV ' . (15) 

i=0 

The corresponding contribution of the periodic orbit pair to the form factor (|^) is then (— l) dv - v '2~ n . 
Note that the amplitudes A v equal 2~ n / 2 for all orbits of period n. 

The from factor can thus be written as a sum over weighted correlations of the degeneracy 
function fbh, i.e. 



K(n,N) = — - 
v ' ; N2 



7, E (^(-^^PnWPnW) 

qeK„(B N ) V q i' I 
qeK„(B N ) V i 7 



(16) 



and [.] denotes the integer part. The sum is taken over the N dimensional integer lattice K n (.B/v) of 
possible vertex staying rate vectors q corresponding to periodic orbits of period n of a binary graph 
-Bjv; the vectors q, q' correspond here to the N/2 components (qi^i, i = 0,N/2 — 1) of the total 
edge staying rate vector only. The contributions of periodic orbit pairs which are not correlated 
by having length differences being a multiple of n are assumed to be Gaussian distributed with 
zero mean. We will neglect these random background contributions from now on and concentrate 
on the contributions from correlated periodic orbit pairs only. 



Before turning to the problem of calculating degeneracy functions, a few remarks on the edge 
staying rates. The components of the edge staying rate vector q are related to each other by the 
shift invariance properties (|l4|). These are N conditions which can be shown to lead to iV — 1 
independent equations for the 2N rates (fo-; together with the restriction 

AT-l 

E « = n ( 17 ) 

i=0 

for orbits of period n, one can write the edge staying rates in terms of N independent quantities, 
which effectively allows to half the dimension of K„(_Bjv). The length degeneracy functions P n 
depends thus on N independent variables only. 

There are further restrictions on the independent components of q. Apart from the obvious 
condition qij > Vij G E(Bn), one must also ensure that the sum over the N independent 
components of q does not exceed n and that the staying rates do correspond to a connected, closed 
path on the graph. An example for an edge staying rate vector q violating the last restriction is q = 
(<Zoo>0, ■ • ■ ,0,Qn-i,n-i) with goo 7^ an d <1n-i,n-i 7^ which corresponds to two disconnected 
periodic orbits. I will come back to the problem of determining the lattice K„ in more detail in 
the next section. 



III. PERIODIC ORBIT LENGTH DEGENERACY FUNCTIONS - ANALYTIC 

RESULTS 

The periodic orbit length correlations in binary graphs with constant transition amplitudes can 
be completely described in terms of the degeneracy function (||). The problem of calculating the 
form factor is thus converted to a combinatorial problem of finding the number of closed (connected) 
paths on a graph which visit each edge the same number of times. This problem can be treated 
explicitly for low-dimensional graphs; results for binary graphs up to order 6 will be presented 
here. 



A. Binary graphs of order N = 2 

The case N — 2 has already been treated by Schanz and Smilansky (1999) in somewhat different 
circumstances. Q We will re-derive some of the results in order to introduce the basic notations and 
concepts which will be useful when considering binary graphs for N > 2. Some new asymptotic 
results for the two-dimensional case will also be presented here. 




FIG. 2. Binary graph of order 2 together with its adjacency matrix A. 

A binary graph of order 2 is shown in Fig. |^. The shift invariance property, Eq. (jlj) , implies 
the following conditions for the edge staying rate vector q = (qoo, goi, Qw, qn), ke. 

go = 9oo + 9oi=9io + 9oo, (18) 
9~i = 9n + 9io = 9oi + 9n 

and go, 9i represent the vertex staying rates. After choosing qoo and qn as independent variables 
and together with the condition ([jj]), one obtains 



1 Schanz and Smilansky (1999) analysed unitary 2x2 matrices in connection with simple quantum (star-) 
graphs. The unitary transfer matrices considered have the extra constraint Lai = Lio. It can, however, be 
shown that this conditions does not lead to additional periodic orbit length correlations, see also Sec. uM. 



qoi = qio = \( n - loo ~ qu), 9o = ~(n + qoo - qu), q\ = ^(n - qoo + qn) , (19) 
for orbits of period n. 

The periodic orbit length degeneracy function P n (qoo, qn) can be derived by starting with the 
special case qoo = qn = 0. One immediately obtains P„(0, 0) = 2 for n even; the two periodic 
orbits correspond to the §— th repetition of the primitive periodic orbits 01 and 10 of period 2. It 
is advantageous to switch to an edge symbol code, i.e., to identify 

00 -> e ; 01 l e 10 -» 2 e ; 11 -> 3 e , 

see also Eq. (|l^ and Fig. ^. The two orbits 01 and 10 can then be written as 

Ie2 e le2 e ■ ■ ■ le2 ej and 2 e l e 2 e l e . . . 2 e l e . (20) 

s / s / 

n n 

The symbol e can only occur after the symbol 2 e and it can be repeated. A periodic orbit of 
period n + qoo can thus be obtained by inserting goo symbols e in between the 2 e l2 blocks in the 
periodic orbit sequences (p0|). Symbols e can be placed at j + 1 positions for the first orbit in ( |2C| ) 
and ^ positions for the second orbit. Similar arguments apply for inserting qu symbols 3 e into 
the sequences ( p0[ ) . Using standard combinatorial formulae to find the number of combinations to 
distribute qoo items among ^ + 1 or | boxes with repetitions, one obtains 

P f„ „ s ( f + 9oo \( f + 9ii - 1 \ , ( f + goo - 1 \f § + 911 

P n+g<m+qil (qoo,qn)=[\ m ){ 2 qu ) + { 2 q o ) { \ xl 

After rescaling (n + qoo + qu) to n and using the relations (|l^) one may write the degeneracy 
function as 

^oo, gil )=^fv i Kv 1 ) • ^ 

9oi V *o /V 9n ) 

The possible integer values for qoo and qu have to obey certain restrictions which follow directly 
from Jig|), i.e. 

9oo + 9n <n and {n - q o - 9n) mod 2 = . (22) 



The degeneracy function fl2l|) approaches a Gaussian distribution in the limit n — > oo; its form 
can be derived with the help of large deviation techniques (Dembo and Zeitouni (1993)), i.e., one 
obtains 

^(9oo,9ii) ~— 2 n e-« 4 * a+ v 2 ) (23) 
Tin 

• i 1 / n \ 1 / \ 

with x= — 7={qoo + 9n ~ t:)> V = — 7=(9oo - Sii) ■ 
n\J2 2 ny2 

The asymptotic result ( p3| ) is too crude to be useful in a calculation of the form factor directly; 
it does provide, however, some insight into the asymptotic behaviour of the various contributions 
entering the form factor. Especially the contributions of topologically degenerate periodic orbit 
pairs, see Eq. (||), can be estimated to be 



K top (n) - J J d qoodqnP^{q 



2 n 

'oo, qn) = — 

irn 



and one obtains a t = log 2 for the growth rate of the diagonal contributions (g) . Periodic orbit 
pairs being degenerate up to a phase difference mir enter the form factor asymptotically as 

{ i) m r f 2 n 2 

K m (n) - 2n+1 / / dqoodqn P n {qoo, qu)P n {qoo + m, qn + m) = e -4 ^ . 



irn 



The form factor thus consists of an increasing number of exponentially growing terms which differ 
in sign (see also Fig. ||). Only a very delicate balance between these terms ensures the cancelations 
necessary to lead to the asymptotic behaviour linv^oo K(n) = 1. The approximations above are 
indeed not sufficient to preserve the asymptotic limit and give exponentially growing terms for large 
n\ similar arguments might hold for the breakdown of semiclassical approximations to quantum 
form factors, see e.g. Tanner (1999). Note also, that the diagonal terms relevant for the HOdA - 
sum rule do not play a prominent role in the discussion above; they give a linear contribution to 
K top which is already sub-dominant for moderate n values. 

The periodic orbit pair contributions to the form factor can be computed explicitly by summing 
the exact length degeneracy function (21) over the possible edge staying rates obtained from the 
conditions (E2[). ft may be written in compact form in the following way (Schanz and Smilansky 
(1999)) 



K(n) 



2 n+l 




(~-l) qol Pn(qoo,qu) 



901=1 



= 1 



(-1) 



n+l 



2 2/+l 



(24) 



with I = [n/2] and qo, qoi can be expressed in terms of 900,911 using jl9|). It is a remarkable 
fact that the sum can be determined explicitly, a result derived by Schanz and Smilansky (1999) 
using quantum graph techniques. The form factor, Eq. (p4|), is displayed in Fig. |3| together with 
the asymptotic results. K(n) approaches 1 in the large n limit, but is different from the RMT- 
result for 2x2 matrices. The periodic structure can be seen to coincide with the start of a new 
family of degenerate orbits and is thus a remnant of non-perfect cancelation of the various K m (n) 
contributions. Convergence of the correlated periodic orbit pair contributions to the RMT - result 
is observed when increasing the order N of the binary graph as will be shown in the next sections. 
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FIG. 3. The form factor for binary graphs of order 2 (dotted line with +) is shown as function of r = ^; 
the partial sums K m contribute with alternating signs starting at r = m + 1; (the dashed lines correspond 
to Km + 1 for m > 0). The GUE - form factor is also displayed for comparison. 



B. Binary graphs of order N — 4 and N = 6 

The edge and vertex staying rates of periodic orbits of a binary graph of order N — 4, see Fig. ||| 
can be written in terms of 4 independent variables. A possible choice for the edge sta yin g rates is 
qoo> qi2, qi\ and 933. The other edge and vertex rates can be computed by using Eqs. (|l4|), explicit 
formulae are given in appendix |a|. 
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110 
11 
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FIG. 4. Binary graph of order 4 together with its adjacency matrix A. The subscript e denotes the edge 
symbol code. 



The periodic orbit length degeneracy function can be obtained by arguments similar to the one 
described in the last section. The discussion is somewhat technical and is referred to appendix [a|. 
The final result is 



™n(900, 912,921, 933) = ~ 
91 



9i 
9oi 



92 
921 



9o-l 
9oo 



93-1 

933 



(25) 



and qi denotes again the vertex staying rates. The possible entries on the 4-dimensional integer 
q lattice can be stated by conditions similar to those in Eq. (|22]). Periodic orbits which differ in 
length by a multiple of 7r have the same vertex staying rates but may differ in the variables 



«o = 9oo + 921, si = gi2 + q 33 . 



(26) 



The length difference for orbits with identical vertex rates is given by AL — ^(Aso + Asi)7r, see 
Eq. (|l5|). The form factor can be written in terms of degenerate periodic orbit pairs only and one 
obtains 



K(n) = - — 
y ' 4 2" 



n-l 



- 1 --go -9i I 

2 + E E 

90+91 = 1 91 = 1 



s +si<n 



E 

i so = 1 900— 921 1, SI = |912— 933 | 



(-l)[^]p„(q) 




(27) 



The form factor K{t) with r = n/4 obtained from Eq. (|27|) is shown in Fig. |6j. It oscillates 
periodically with decreasing amplitude about the RMT - result similar to the behaviour observed 
in the case N = 2, see Fig. |[ A closed expression for the sum similar to Eq. ( p7j ) could not be 
found. 

The sums in d27j ) are already quite cumbersome and the number of summation variables increases 
with the order N. The number and complexity of the restrictions for the q - lattice JK. n (Bw) in- 
creases accordingly. The case N = 6 can, however, still be treated along the ideas developed 
above; it will be presented here as a last example for obtaining the form factor by summing over 
the periodic orbit length degeneracy function. 
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FIG. 5. Binary graph of order 6 together with its adjacency matrix A. The subscript e denotes the 
symbol code. 

The binary graph of order N = 6 is shown in Fig. ||. A possible choice for the independent edge 
staying rates is goo, 9i3, 924, 931, 942, and 955. The derivation of the degeneracy function can again 
be found in appendix [A], the final result is 



nqw 

9293 



31 

531 



92 
942 



93 

gi3 



94 
924 



9o-l 
9oo 



95-1 

955 



(28) 



The vertex rates q~i and the edge rate q\2 entering (|28|) can be expressed in terms of the independent 
variables q = (goo, 913, 924, 931, 942, 955), see appendix [A]. The summation over the six dimensional 
lattice K„(_Bg) of possible q vectors can be stated in terms of the vertex staying rates and the 
variables 



«o = 9oo + 931, si = 912 + 943, s 2 = 924 + 955 ■ 
The expression for the form factor as sum over degenerate periodic orbit pairs is thus 

n [(«— 9o)/2] [(n— qo— 2gi)/2] f (s +si+s 2 )<« 



^)4iE E 



90 = 1 gi=0 



E 

92=0 



E (-D 



Pn(q) 



(29) 



and the inner sum runs over all possible s,, i = 0, 1, 2 values. The form factor if (r) after summing 
Eq. @ is displayed in Fig. | with r = ra/6. The sums, Eqs. @ and @, follow the RMT - 
result more closely than in the N — 2 case, see Fig. 0. The linear behaviour for r < 1 starts to 
emerge and convergence to the asymptotic result K — > 1 is observed in the large t = n/N limit. 

Larger matrices have to be considered in order to test convergence of degenerate periodic orbit 
pair contributions towards the RMT form factor for all r. Determining the degeneracy function 
and the lattice conditions K„(i3/v) becomes increasingly difficult for graphs of order N > 6. In the 
next section, I will therefore present results obtained from counting all correlated periodic orbit 
pair contributions directly. 
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FIG. 6. Periodic orbit pair contributions to the form factor for binary graphs of order N — 4 and TV = 6. 



IV. PERIODIC ORBIT PAIR CONTRIBUTIONS TO THE FROM FACTOR FOR DE 

BRUIJN GRAPHS OF ORDER N > 8. 

The periodic orbit pair contributions to the form factor can be calculated directly by determining 
the set of periodic orbits of given period n and calculating periodic orbit degeneracies with the 
help of edge and vertex staying rates and the condition ([l5|). The task of finding the set of periodic 



orbits is especially simple for de Bruijn graphs, i.e. for binary graphs of order N = 2 r , due to the 



one-to-one relation between periodic orbits and finite binary symbol strings, see section II B 



Counting the periodic orbit pair degeneracies explicitly does, however, seriously limit the range 
of periods over which periodic orbit correlations can be considered. Due to the exponential increase 
in the number of orbits only values up to n ~ 26 could be reached. This in turn sets an upper 
bound on the r = n/N values for which the form factor can be studied. 



CUE - results : 

Results for N = 8, 16 and 32 and no further symmetry present are shown in Fig. [?]. One observes a 
convergence of the periodic orbit pair contributions to the GUE-result; the kink at r = 1 is resolved 
for binary graphs of order N = 16; the periodic orbit results follows the linear behaviour for r < 1 
even closer for N — 32. It was not possible to extend the results for n — 32 to the critical time 
r = 1 due to the restrictions on the available n values. The small r behaviour is dominated by the 
exponentially increasing topological contributions, see Fig. 0. The so-called diagonal contributions 
due to cyclic permutations of periodic orbit codes are important only in the small t regime, i.e. 
t ~ log 2 (N)/N, before vertex exchange degeneracies set in. 




GOE - results : 

So far only unitary transfer matrices without symmetries have been considered. Symmetries in the 
dynamics impose additional correlations on periodic orbit length spectra and do have an effect on 
the spectral statistics. Time reversal symmetry is of special importance as it occurs frequently in 
physical systems; correlations due to time reversal symmetry are in addition non-trivial leading to 
a form factor which is not piecewise linear as in the GUE case; only the linear behaviour for K(t) 
in the limits r — > and r — > oo is understood in terms of semiclassical arguments (Berry (1985)). 

It is a priori not clear how to establish time reversal symmetry for the dynamics on an arbitrary 
directed graph. Time reversal symmetry can, however, be constructed for de Bruijn graphs of order 
N = 2 k using the underlying binary symbolic dynamics and the edge code, Eq. (|l0|). The edge 
code can be written in terms of a binary string of length k + 1 such that i e — J2i=i a i(*e)2 fc+1 ~', 
and a(i e ) = (ax, . . . <2fe+i) is a binary string of length k + 1 with a; G 0, 1. Time reversal symmetry 
can be established by imposing 

L ie =Li* e , r ic = r,v if a(i e ) = a.(i' e ) (30) 



for the edge lengths and transition rates and a denotes the code a written backwards, i.e. 



'}2d=i a ((*e)2 i_1 - The condition, Eq. (30), and the unitarity condition, Eq. (|l2|), are assumed to 
be the only sources of correlations in the periodic length spectrum. 

Time reversal symmetry does not effect graphs of the order N < 8, N = 2 k . This is due to 
the fact that the edge staying rates for a given edge and its time reversed partner are related by 
conservation laws ([14|) such that there are no further degeneracies for these low dimensional cases. 
One finds for N = 2, for example, that goi = qio and a periodic orbit and its time reversed partner 
are always in the same degenaracy class even if the condition ( |30| ) is not imposed. 



K 




FIG. 8. Periodic orbit pair contributions to the form factor for binary graphs of order N — 8, N = 16 
and N = 32 and time reversal symmetry. 



Results for graphs with time reversal symmetry are shown in Fig. pr the case N = 8 is indeed 
identical to the non-time reversal symmetric result in Fig. 0. The results for N — 16 and 32 
are, however, different from those in Fig [?]; the periodic orbit pair contributions approach the 
GOE - result and not the GUE form factor with increasing TV. The condition ( po|) does therefore 
introduces new correlations among periodic orbits for N > 8 which are beyond the additional 
topological degeneracy between an orbit and its time reversed partner giving rise to a factor 2 in 
the HOdA - diagonal approximation. Note also the exponentially increasing components for small 
t due to topological degeneracies similar to those in Fig. 0. 



V. CONCLUSIONS 



Degeneracies in the length spectrum of periodic orbits of generic directed graphs have been 
studied. Transition rates and edge lengths in the graph are identified with the amplitudes and 
phases of matrix elements of the complex transition matrices. General concepts like edge and 
vertex staying rates as well as the periodic orbit length degeneracy function have been introduced. 
The form factor can be written in terms of the degeneracy function revealing an exponentially 
increasing 'diagonal contribution' due to topologically degenerated orbits. Topological degeneracies 



2 The unitary 2x2 matrices studied by Schanz and Smilansky (1999) do thus correspond to time reversal 
symmetric (binary) graphs. 



exist independently of the actual choice of length segments on the graph (defined through the 
transition matrix) and are a purely 'classical' effect depending only on the topology of the graph. 
Further correlations amongst orbits are introduced when considering unitary transfer matrices. 

These correlations have been studied for a particular simple class of graphs, so-called binary 
graphs with constant transition amplitudes. The correlations can be given explicitly in terms 
of edge and vertex staying rates. One finds in particular that periodic orbits which have the 
same vertex staying rates differ in length by exactly a multiple of tt. Finding the periodic orbit 
degeneracy function turns into a combinatorial problem which has been solved for binary graphs 
with up to 6 vertices. 

The form factor can be shown to consist of exponentially increasing contributions which balance 
in a very delicate way to give lim^oo K(t) — 1. The periodic orbit sums also reveal convergence 
towards the RMT - result for intermediate r - values when increasing the order of the graph, both 
for time reversal and non-time reversal symmetric binary graphs. Binary graphs may thus turn 
out to be an ideal model systems to study the connection between periodic orbit formulae and 
random matrix theory. Eigen spectra of binary graphs follow generic random matrix statistics in 
the large N limit and periodic orbit correlations are known explicitly 
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APPENDIX A: PERIODIC ORBIT LENGTH DEGENERACY FUNCTION FOR 
BINARY GRAPHS: EXACT RESULTS 

The expressions for periodic orbit length degeneracy functions for binary graphs of order N = A 
and N = 6, Eqs. ( p5[ ) and (|28|), will be derived here. 

The case N = 4: 



As for binary graphs of order N — 2 discussed in section III A, it is useful to switch to an edge 
symbol code, see Eq. (|To|); adopting the vertex symbol code of Fig. |J, one defines the edges as 

00 ^0 e ; 01->l e 12 ^2 e ; 13 -> 3 e ; 
20 -► 4 e ; 21 -► 5 e 32 -► 6 e ; 33 -> 7 e . 

A suitable set of independent edge staying rates is qo c , q2 e , 95 e , 97,. and I will drop the subscript e 
as long as there is no confusion with the vertex code. The remaining edge and vertex staying rates 
for periodic orbits of period n can be written in terms of the edge staying rates above with the 
help of Eqs. (|l4|), i.e., one obtains 

9i = 94 = 7 (n - qo + 92 - 3q 5 - q 7 ) (Al) 
93 = <?6 = T (n - qo - 3q 2 + 95 - qi) 

for the edge rates and 

qo = ^ (n + 3q + q% - 3^5 - qi) (A2) 

1 . 

<7i = qi = ^ (n - go + 92 + 95 - 97) 

93 = ^ (n - 9o - 3<?2 + 95 + 3q 7 ) 

for the vertex staying rates qi and the index i denotes the vertex code, here. 

The periodic orbit length degeneracy function P n (qo, 92, 95, 97) can be computed by starting from 

n n 

-Pn(°7 2 1 ^>o) = 2 f° r n even ; 

the set of edge staying rates above corresponds to the orbits of length n with edge symbol code 

2525. ..25 and 525. ..252. (A3) 

One proceeds by noting that an edge symbol '2' in the sequences (|a|) can be replaced by the 
sequence '3 6' to give an orbit of length n + 1. Similarly one may substitute a symbol '5' by the 
sequence '4 1'. Replacing m symbols '2' and k symbols '5', m,k < one obtains 

w fc (o,f-m,|-^o)=2(|)(f) + (|: i 1 )(|) + (n(f: 1 i 



and the last two terms in the sum come from orbits which start with a symbol '6' or a symbol '1', 
respectively. After replacing n + m + k by the new periodic orbit length n', i.e., n = n' — m — k 
and writing qi = j(n' + qi + q§) with qi = \(n' — 3m — k), q$ = \ (n' — m — 3k) one obtains 



p„.(o, 92 , 55 ,o) - 2 ( * ) ( % ) + ( * : ; ) ( % ) + ( 1 ) ( %-_\ ) . (A4> 

Next, one notes that an edge symbol '0' or '7' can be inserted between any symbol '4' and '1' or 
'3' and '6', respectively, to obtain a periodic orbit of length n' + 1. Inserting qo symbols '0' and q 7 
symbols '7' into fc sequences '4 1' and m sequences '3 6' (with repetition) leads to 



Pn'+ qo+q7 (qo,q2,q5,q7) = 2 ( ^\ 



qi\ ( qi\ ( m + q 7 -l \ ( k + qo-1 

qi J \ qo 

qi - 1 \ / qi \ ( m + q 7 \ ( k + q - 1 
m-1 J \ k J y q 7 J \ q 

q%\ ( q% - 1 A ( m + q 7 -l \ ( k + q a 
m J \ k-l J \ q 7 J \ q 

After rescaling to the new periodic orbit length n" = n' + qo + q 7 and summing the three contri- 
butions, one obtains 

^•.«).^(*)(t)(" + :-)(* + r I ) m 



with qi = \{n" — qo + qi + 95 — qi) as in (A2). The final result (|25| ) is obtained after noting that 
m = <73 = qo = q\ — qi and k = q\ = q^ = q% — q$. Special care has to be taken in the case qo or 
97 = 0. 

The case N = 6 : 

The periodic orbit length degeneracy function for N = 6 can be derived by ideas similar to those 
outlined for N = 4; I will sketch the main steps here and leave the details to the reader. 
An edge symbol code is defined starting from the vertex symbol code used in Fig. || to be 



00 - 


+ e ; 01- 


■* le 


12- 


+ 2 e ; 


13- 


y 3 e ; 


24- 


^4 e ; 25- 


->5 e 


30 - 


■+ 6 e ; 


31 - 




42 - 


^8 e ; 43- 


^9 e 


54- 


■+ 10 e ; 


55 - 


-> He 



A suitable set of independent edge staying rates is qo s , 93 e , 94 e , 97 c , 9s c , 9n e and I will drop the 
subscript e from now on. The other edge staying rates of orbits of period n are then given by 

qi = 96 = ~ (n - q + 3q 3 + q 4 - 5q 7 - 3q 8 - q n ) (A6) 
6 

qi = g (n - 9o - 3g 3 + 94 + 97 - 3q 8 - gu) 
95 = 9io = g (n - 9o - 3<? 3 - 5g 4 + 97 + 3<? 8 - 9n) ; 

the vertex rates are 

9o = 9o+9i; 91 = 93 = 92 + 93 ; 92 = 94 = 94 + 955 95 = 9io + 9ii- (A7) 

A suitable starting point for the periodic orbit length degeneracy function is the periodic orbit 
'2 4 9 7' (in edge code) or '1 2 4 3' in vertex code, see Fig. g. One obtains 

Ti TL 

P„(0,0,g 4 = -,97 = ->°>°) = 4 for nmod4 = 0. 



A symbol '7' can be followed by a loop '3 7' (with repetition), a symbol '4' may be followed by a 
loop '8 4' (with repetitions). Inserting k loops '3 7' and m loops '8 4' into a sequence '2 4 9 7...' 
of length n — 2k — 2m yields 



P n (0, k, hji - 2fc + 2m), i(n + 2fc - 2m), Z, 0) = 

,, . \{n + 2k - 2m) - 1 \ / \{n - 2k - 2m) - 1 \ / ±(n + 2fc - 2m) \ ( \{n - 2k - 2m) - 1 
fc y I m J \ k J \ m 

+ I |(n + 2fc - 2m) - 1 \ / \(n - 2k - 2m) \ + / |(n + 2fc - 2m) - 1 ^ ^ i(n - 2fc - 2m) - 1 



m 



+ 



i(n + 2fc - 2m) - 1 \ / \{n - 2k - 2m) - 1 
k ) [ m- 1 



and the different terms in the sum correspond to a first symbol in the periodic orbit code being 
either '2' or '9', '7', '4', '3' or '8', respectively. Next, one notes that every symbol '7' or '4' can 
be replaced by the sequence '6 1' or '5 10', respectively. I omit the somewhat lengthy combinato- 
rial expressions here. The full periodic orbit length degeneracy function is finally obtained after 
inserting symbols '0' or '11' into the sequences '6 1' or '5 10', respectively, and summing over the 
various binomial coefficients. 



